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ABSTRACT 

Recent studies have shown that the local mass-metallicity relation depends on the specific star forma¬ 

tion rate (SSER). Whether such a dependence exists at higher redshifts, and whether the resulting M*-Z-SER 
relation is redshift invariant, is debated. We re-examine these issues by applying the non-parametric techniques 
of Salim et al. (2014) to ~ 130 z ~ 2.3 galaxies with N2 and 03 measurements from KBSS (Steidel et al. 2014). 
We find that the KBSS M*-Z relation depends on SSER at intermediate masses, where such dependence exists 
locally. KBSS and SDSS galaxies of the same mass and SSER (“local analogs”) are similarly offset in the BPT 
diagram relative to the bulk of local star-forming galaxies, and thus we posit that metallicities can be compared 
self-consistently at different redshifts as long as the masses and SSERs of the galaxies are similar. We find that 
the M*-Z-SER relation of z ~ 2 galaxies is consistent with the local one at logM* < 10, but is offset up to 
-0.25 dex at higher masses, so it is altogether not redshift invariant. This high-mass offset could arise from 
a bias that high-redshift spectroscopic surveys have against high-metallicity galaxies, but additional evidence 
disfavors this possibility. We identify three causes for the reported discrepancy between N2 and 03N2 metal¬ 
licities at z ~ 2: (1) a smaller offset that is also present for SDSS galaxies, which we remove with new N2 
calibration, (2) a genuine offset due to differing ISM condition, which is also present in local analogs, (3) an 
additional offset due to unrecognized AGN contamination. 

Keywords: galaxies: evolution—galaxies: fundamental parameters 


1. INTRODUCTION 

The relationship between the stellar masses, gas-phase 
metallicities, and star formation rates of galaxies (M*-Z- 
SER) has received increasing attention over the past ^5 years, 
both because of its conceptual simplicity and its potential to 
provide deep insight into the processes that regulate galac¬ 
tic star formation (SE) and drive galaxy evolution over cos¬ 
mic time. While the local (z < 0.3) luminosity-metallicity 
and mass-metallicity ( M„-Z) relations have b een studied for 
decades, beginning with|Lequeux et al.|(|1979|l, only more re¬ 
cently has the star formation rate (SLR) been proposed as a 
second parameter in the local mass-metallicity relation ( |E1-| 
lison et al.||2008|). The primary^ correlation shows metallic- 
ity increasing with stellar masj^unt il a plateau is reached at 
log M* ~10.5 (|Tremonti et al.|2004|l, while the general sense 
of the secondary dependence with SER is that at fixed mass, 
galaxies with higher SERs tend to be more metal poor ([Elli¬ 


son et al.|2008[|Mannucci et al.|2010]|Lara-L6pez et al.|2010[ 

Hunt et al.|2012|l. ' ' 


Conflicting results on the characterization of the local M*- 
Z-SER relation have ignited debate over whether the reported 
secondary dependence o n SER could be spurious, and due to: 
sample selection effects (|de los Reyes et al.|20151l, correlated 
errors in the measurements of SLR and metallicity (Lilly et al. 
2013|l, systematic errors in metallicities (|Yates et al.|2012|l, or 
biases introduced by SDSS fiber spectroscopy (|Sanchez et al.' 
2013|l. In|de los Reyes et al.|(|2015|l, the limitations of previous 
parameterizations used to characterize the local relation and 
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the need for non-parametric techniques to enable comparative 
analysis of different datasets were also highlighted. 

Thus, to gain insigh t into t he origins of the conflicting re¬ 
sults, in |Salim et al.| (|20T4|l we devised a non-parametric 
analysis framework based on the SER offset from the local 
star-forming (“main”) sequence at a given M*, and under¬ 
took a comprehensive re-analysis of the local M»-Z-SER re¬ 
lation using a Sloan Digital Sky Survey (SDSS) dataset to¬ 
gether with GALEX ultraviolet and WISE infrared photome¬ 
try. Studying the M*-Z-SER relation in terms of SERs or 
specific SERs (SSERs) relative to typical values on the “main 
sequence” is more physically motivated t han using absolute 
SERs , which, to first order, scale with M* ( |Brinchmann et al.] 
|2004| l. Although we concluded that the dependence on SSER 
IS not spurious (after investigating multiple SER and metal¬ 
licity indicators), the analysis e xposed important fe atures of 
the relationship. In particular, |Salim et al.| P014)l showed 
that adding the SER as a second parameter does not greatly 
decrease the scatter in the M*-Z relation when the metallici¬ 
ties of individual galaxies, rather than the median-binned val¬ 
ues are considered (i.e., the M*-Z-SER relation is not tight). 
We confirmed that the overall SER dependence is weaker, or 
absent, at higher masses. However, at a given mass the de¬ 
pendence on SER is much stronger for intensely star-forming 
galaxies above the “main sequence” (galaxies with high SSER 
for their mass). We noted that simple parameterizations of the 
local M*-Z-SER relation (a plane, or the projection of least 
scatter in Z) do not capture this behavior of galaxies with high 
relative SSER because the parametrizations are dominated by 
the bulk of “normal" galaxies along the core of the main se¬ 
quence. 

Recognizing the limitations of parameterizations of the lo¬ 
cal M*-Z-SER is particularly important in the context of test¬ 
ing its redshift invariance. The concept of invariant M ^-Z- 
SFR relation was put forward by|Mannucci et al.|(|20101l, who 
refer to a it as the “fundamental metallicity relation” (FMR). 
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They showed that selected galaxy samples up to z ^ 2.2 lie 
along the projection of the local M*-Z-SFR that minimizes 
the scatter in metallicity, i.e., that they are consistent with a 
non-evolving M*-Z-SFR relation. M*-Z relations at z > 0 
then represent the slices of the invariant (i.e., fundamental) 
M*-Z-SFR relation at (S)SFRs higher than those found in lo¬ 
cal galaxies. However, a number of more recent studies have 
concluded that samples do not lie on the local 
SFR relation, implying an evolving relation. M ost of these 
studies were based on |Mannucci et ar](2010|l parametriza- 
tion of the local relation (e .g., [Cullen et al.|2014[|Zahid et al.| 
2014bt |Maier et al.||2014|l, except for|Sanders et al.| (|2015|l, 
who perform ed a airect compariso n \vith the local sampled 
The result of|Sanders et al.|(|2015|l highlights other possible 
causes for the discrepant results, e.g., an evolution in metallic- 

Steidel et al.|20T4ll. 



may also play 

Establishing the existence of SFR dependence of M*- 
Z relation at higher redshifts would be important even if the 
resulting high-redshift M*-Z-SFR relation did not coincide 
with the one followed by local galaxies. The evidence that 


elusive (jCresci et al. [2012 

[Zahid et al.j[2014bl [Wuyts et al. 

2014 

Steidel etal.|2014 M 

laier et al. 2014 de los Reyes et al. 


Sanders et al.|2015[l. 


Armed with a more detailed, non-parametric characteri¬ 
zation of the local M*-Z-SFR relationship we can examine 
these points of contention in recent work involving galaxies 
at z 2. Here, we will apply our analysis framework to re¬ 
examine whether a secondary dependence of the M*-Z rela¬ 
tion on SFR also exists for higher redshift samples, and more 
generally, whether the local M*-Z-SFR relation is invariant 
with redshift and describes star-forming galaxies at all stages 
of their evolution over cosmic time. We focus on the latest 


and most comprehensive sp ectroscopic datasets (Steidel 
et al. 112014} Sanders et al. 2015|l. We present evidence that 


this approach is able to circumvent the thorny issues of incon¬ 
sistent or evolving metallicity calibrations. We highlight the 
possib le role of observational selection effects ( ] Juneau et al.| 
|2014| l and/or the possibility that some high-mass high-redshift 
galaxies contain unrecognized contribution from AGN line 


emission. 


2. DATA AND SAMPLES 


Our analysis is p rimarily based on the dataset published in 
Steidel et al.|P014) l from the Keck Baryonic Structure Survey 
a near-lR spectroscopic survey performed with the 
MOSFIRE multi-object slit spectrograph on the Keck I tele¬ 
scope. KBSS has obtained H and K-band spectra of galax¬ 
ies at 1.95 < z < 2.65. The majority of targets had previ¬ 
ously determined redshifts from optical spectroscopy and had 
been originally selected u sing a variety of techniques based 
on rest-fr ame UV colors (|Adelberger et ^ 2004[ [Steidel et| 


l.|(|2014|l analysis includes indi- 
TNTatlos of [OIII]5007 and Ha 


ai:i[2004ll. The [Steidel et al.l 
vidual galaxies where the S/f 
line fluxes exceed 5, and those of [NII]6584 and H/3 exceed 
2. For an average exposure time of 11000 s, these cuts cor¬ 
responds to Ha and [OIII] 5cr limits of ^ 10“*' ’ e rg s~* (ob¬ 
tained from scaling MOSDEF s ensitivity given by [Coil et al.[ 
[2015J. The published version of[Steidel et al.[([20l4p presents 
161 galaxies with both N2 = lc)g[F([NllJ6584)/F(Ha)] and 
03 = log[F([Oni]5007)/F(H/3)] line ratios and an additional 
31 galaxies with only the N2 measurements (i.e., currently 
lacking //-band observations). The total of 192 galaxies ex¬ 


cludes seven objects identified as AGN based on broad emis¬ 
sion lines or the presence of higher ionization species. Stellar 
masses were derived from SED fitting, and SFRs from Ha 
fluxes, corrected for slit losses and for extinction based on the 
continuum dust attenuation estimate from the SED fits. We 
note t hat the data used in o ur analysis, however, are the subset 
of the[Steidel et al.[([2014|l sample, which was presented in the 
preprint version of their paper. We use the smaller, preprint 
sample because the published version of the paper omitted 
SFRs from the tables. Thus, the dataset used here contains 
108 galaxies with both N2 and 03, and an additional 18 with 
N2, i.e., 2/3 of the published sample. Line ratios and stellar 
masses are taken from the published tables, and SFRs from 
the preprint. Differences between preprint and published line 
ratios are not significant (0.07 dex scatter, with no system¬ 
atic offsets), so we assume that the SFRs have not changed 
significantly either, as also evidenced by a similar appearance 
of figures that involve SFRs in the preprint and the published 
paper. 

We augment the analysis with measurements from the 
MOSFIRE Deep Evolution Field (MOSDEF) survey early ob- 

et al.| 


servations ( Kriek et al.|2014[[Coil et al.|20151 [Shapley et 
2015| l. MOSDEF targets //-band selected galaxies at 2.1 


z < 2.6 within the CANDELS survey areas. Targeting priori¬ 
tization criteria include the availability of a spectroscopic red¬ 
shift from previous work (40% of the sample), brightness, and 
photometric redshift in the target range. 03N2 (= 03-N2) 
and N2-based metallicities were obtained for 53 galaxies for 
which all four lines have S/N ratio > 3, c orresponding to H a 
and [OIII] 3cr limits of ~ lO'**'^ erg s"^ ( jCoil et al.poisf 
AGNs were removed from t his sample based on X-ray and IR 
indicators, or if N2> -0.3 ( [Coil et al.|2015| l. Stellar masses 
come from SED fitting, and SFRs from Ha fluxes corrected 
for extinction using the Balmer decrement. Data tables with 
measurements for individual MOSDEF galaxies are no t yet 
published. H owe ver, we used the publi shed figures in Coil| 
[et al.[ ([2015|l and [Sanders et al.[ ([2015J to compare M»-Z, 
SSFR-M* and 03-N2 relations from MOSDEF with those 
from SDSS and KBSS. 

Comparison with the relations followed by loc al galaxies 
is based on the SDSS DR7 spec troscopic sample ( [Strauss ^ 
[^[2002[ Abazajian et al.|2009[ l processed by the MPA/JHU 
group. Sample selection follows that of [Mannucci et al.| 
([2010[l, except that we extend the redshift range (0.005 < 
z < o! 3) as long as the mass included in the spec troscopic 
fibe r is > 10%, following what we have done in Salim et| 
[ak] ( [2014| l. Dropping the redshift limit from 0.07 to 0.005 
removes low-SFR inco mpleteness for star- forming galaxies 
with 9 < logM* <10 ([Salim et al.[[2014[). Typical mass¬ 
covering fraction of SDSS fiber spectroscopy is 30%, with 95 
percentile range between 17% and 50%. Inclusion of lower- 
redshift galaxies, which have lower mass-covering fraction, 
does not affect any of the results. It is also important to note 
that the selection is based solely on S/N ratio in Ha being 
above 25. The limit on only the Ha line ensures that the sam- 
ple is not biased in metallicity ([Salim et al.[[2014| [Mannucci] 
et al.|2010|i , while it is high enough that other required nebu¬ 


lar lines will be well measured. Lowering the limit to values 
below 25 increases the number of galaxies with lower rela¬ 
tive SSFR. However, these galaxies entirely follow the trends 
established by galaxies with higher S/N ratio, just with less 

^ While KBSS is on average deeper than MOSDEF, the former requires 
higher detection threshold in [OIII] line. 





























































































































MASS-METALLICITY-SFR RELATION 


3 




Figure 1. The impact of the choice of metallicity indicator on inferred evolution in the M*-Z relation. Local (z ~ 0) sample comes from SDSS (green hand 
is the 95 percentile range at a given mass; white line shows an average trend), while z ~ 2 data are from KBSS (black points; binned averages shown as red 
lines) and MOSDEF (binned averages shown as blue lines). MOSDEF averages are based on galaxies with M* > 10, for which their sample is not affected by 
incompleteness. Different average lines are obtained by v arying the starting posi tion for binning in fifths of the 0.25 dex bin width. Left (middle) panel is based 
on N2 (03N2) line ratios, converted into metallicity using|Pettini & PagelH2004I linear calibrations. Right panel shows N2 after the recalibration is performed in 
order to bring N2 and 03N2 metallicities into a better mutual agreement tor the SDSS sample. 03N2 shows a higher degree of “evolution” than N2, especially 
at higher masses. MOSDEF and KBSS MZRs agree between each other for N2, but depart to some extent for the more massive galaxies in the case of 03N2. 


precise metallicity measurements. 

For init ial analysis, we exclude SDSS galaxies that lie 
above the [Kauffmann et al. ( |2003[ ) AGN demarcation line. 
We use total (integrated) SFRs and stellar masses from the 
MPA/JHU catalog , wh ich are determined following |Brinch-| 
|mann et ah] ( |20()4 1 and |Salim et al.|p007| ), respectively, with 
additi on Jdetails given i n onli ne documentationp] Specifi¬ 
cally, |Brinchmann et al.| (|2004|l SFRs are based on a hybrid 
combination of an emission-line based SFR within the spec¬ 
tral fiber and a photometric estimate outside of the fiber, and 
hence should capture the total activity of the galaxies. 

A Chabrier IMF and standard cosmology (Hq = 70 km s“' 
Mpc“', ilm = 0.3, r^A = 0.7) is assumed throughout the paper. 

3. METHOD 


As in Salim et al. ( 2014) l, we examine the metallicity as a 
function of the relative SSFR at fixed stellar mass. We define 
the relative SSFR as the offset from the local (SDSS) star¬ 
forming (“main”) sequence; 


dependence on the mass and SFR to be simult aneously ex- 
plored, in contrast to a plane parameterization of|Lara-L6pez| 
et al. (|2010|l, which forces a fixed dependence on both the 


mass and the SFR, or the projection of least scatter in Z of 
[Mannucci et al.| ( |2010| l, which forces a fixed SFR dependence 
at a given mass. Moreover, our methodology allows one to 
test whether high-redshift samples follow the same M*-Z- 
SFR relation as the local galaxies without extrapolation of 
the assumed parametrization into regions which may not be 
well populated by local galaxies, and therefore do not carry 
much weight in the parameterization, but are occupied by 
high-redshift galaxies. 

We will also apply a related non-parametric method for an 
even more direct test of M*-Z-SFR invariance. The method 
consists in comparing the metallicities of z ^ 2.3 galaxies with 
the metallicities of their local “analogs”. A local “analog” is 
defined as an SDSS galaxy with and SSFR most similar 
to a given high-redshift galaxy, i.e., a galaxy for which the 
metric: 


AlogSSFR = log SSFR- (logSSFR)M., (1) 

where (logSSFR)M. is the median log SSFR of galaxies with 
M*. Alternatively (not used here), (logSSFR) m, can be a 
value obtained from fitting the local “main” sequence with 
some mean relation (e.g., 

(logSSFR)M. =-0.35(logM*-10)-9.83 (2) 

from [Salim et al.|[2007| l. We define relative SSFRs in refer¬ 
ence to local main sequence (as opposed to high-redshift one) 
because the local SSFR-M* relation is more robustly known. 
In principle, relative SSFRs can be defined with respect to 
high-redshift sequence as well, and this change in zero point 
would not impact the analysis. 

Using this simple analysis framework, one can deter¬ 
mine if an SFR dependence is present without assuming 
a parametrization of the M»-Z-SFR relation. Such non- 
parametric techniques are essential because they allow the 


http://www.mpa-garching.mpg.de/SDSS/DR7 


D^ = {AlogSSFRf + (AlogM^f (3) 

is minimized. Large volume of SDSS allows finding very 
close matches in SSFR and M* (small D). We use the word 
“analog” with caution, because we are not implying that such 
local galaxies undergo the same physical processes as the 
high-redshift galaxies, but only that in the context of invariant 
M*-Z-SFR relation, the metallicities of high-redshift galax¬ 
ies and the local “analogs” should be the same. 


4. MAIN RESULTS 

4.1. Mass-metallicity relations and the systematic 
differences in local metallicity calibrations 

To set the stage for subsequent analysis, we start by in¬ 
specting the Mtf-Z plots (Figure [T]l, which compare the lo¬ 
cus of SDSS galaxies with those from KBSS and MOSDEF 
samples. Left and middle panels of Figure [1] are based on 
N2 and 03N2 line ratios conver t ed into metamcity using the 
calibrations of jPettini & Pagelj (|2004|), which were derived 
using mostly direct-m ethod abundance s of HI I regions in lo¬ 
cal galaxies. For N2, jPettini & Pagelj (j2004j) provide linear 
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Figur e 2. Comparison of metallicities of SDSS galaxies derived based on 03N2 and N2 line ratios. 03N2 metallicity in all panels is derived from|Pettini &| 
|Pagel|ij20(M^ calibration, while N2 metallicity is based on lineai* and cubic|Pettini & Page lH2Q04) relations (left and middle). One expects IPettini & PagelH2UU4f 
calibrations, being local, to yield consistent metallicities for SDSS galaxies, but this is not entirely the case. We therefore recalibrate N2 calibration by assuming 
a rational functional form (Equation!^ and minimizing the deviations with respect to 03N2 metallicities. The resulting new N2 metallicities show no systematic 
offsets (right panel) and we use themior the remainder of the paper. 


and cubic calibrations. Figure [T] (left panel) us es linear cal¬ 
ibration following KBSS and MOSDEF studies( |Steidel et aE] 
|2014[ [Sanders etal.|2015[ ). 

The first point to be emphasized from Figure is that the 
M*-Z relations based on N2 and 03N2 metallicity indicators 


are not consiste nt even for SDSS galaxies, even though Pettini 
|& Pagel| ( |2004| ) calibrations have been based on local galaxies. 
Most notably, N2 metallicities at higher mass (logM* > 10) 
are 0.11 dex lower than 03N2 metallicities. There are several 
reasons for the local mismatch, which can be best appreciated 
from Figure® where the metallicities of SDSS galaxies based 
on 03N2 and N2 indicators are compared directly. First, the 
linear N2 calibration is too crude to capture the saturation of 
N2 at high metallicities, and thus leads to a diverging offset 
(Figure]^ left). This point is relevant for 12+log(0/H)> 8.6 
and thus mostly affects the local samples which contain galax¬ 
ies with such super-solar metallicities. The relationship be¬ 
tween N2 and metallicity in this regime is somewhat better 
describ ed with the cubic N2 calibration of [Pettini & Page!] 
( [2004| ), but there are still signifi cant systematics at high m etal- 
hcities (Figure [^ middle; also [Kewley & Ellison|2008'[ their 
Figure 2). Second, [Pettini & Pagel ( ®)04p calibrator saniple is 
relatively small (10"^), so some differences due to the limited 
accuracy of the functional fitting are to be expected when the 
comparison is performed within a much larger ( 10^) SDSS 
sample. Furthermore, the [Pettini & Pagel[ (|2004|l calibrator 
sample may not be entirely representative of a typical low- 
redshift population. Consequently, even for lower metallici¬ 
ties (12H-log(0/H)~ 8.4), the metallicity based on linear cal¬ 
ibration of N2 is somewhat offset with respect to 03N2—it 
is ^ 0.03 dex higher (Figure [^ left). An offset of a simi¬ 
lar magnitude, but in the opposite direction, is present for the 
cubic N2 calibration (Figure [^ middle). This discrepancy at 
lower metallicities will be important for high-redshift galaxies 
where measured metallicities tend to be lower. Thus, the 0.13 
dex systematic diffe rence between N2 and 03N2 metallicities 
of z ~ 2.3 samples ( [Newman et al.|20T4 [Zahid et al.|2014b[ 
Steidel et al.|2014|, which is usually attributed to evolution in 


one or both calibrations, is actually in part (~ 1 /4) due to the 
mismatch in local calibrations. 

In order to be able to separate possible systematics arising 
from evolution (i.e., the inapplicability of local calibrations at 


high redshift) from those that stem merely from mismatch of 
the local metallicity scales, we recali brate the N2 relation to 
match the metallicities resulting from [Pettini & Pagel| (|2004|l 
03N2 calibration. We take 03N2 metallicities as hducial be- 
cause the 03N2 calibration is not subject to saturation, so it is 


more likely to b e close to the linear form assumed in Pettini 
[& Pagel[ ( [2004[ ). The new N2 calibration is represented with 
a rational function that tends to a linear relation for low val¬ 
ues of N2 (low metallicities), and has asymptotic behavior for 
high values. Its analytical form is: 


12-l-log(0/H)N2 = a+^N2-l-c/(N2-t/) (4) 

and max(12-l-log(0/H)N2) = e, where a,b,c,d and e are the 
free parameters obtained from the minimization of the sum of 
the binned median deviations of orthogonal offsets of points 
in Figure [^from the diagonal (the 1:1 relation). Minimization 
is performed on binned values in order to give uniform weight 
at a range of metallicities. Parameter d represents the posi¬ 
tion of the asymptote in N2, while e is the value above which 
metallicities are assigned the value e. Parameter e signifies 
the point above which N2 cannot provide reliable information 
on metallicity (except that it is high). The relation with the 
best-fitting parameters is: 


12-i-log(0/H)N2 = 8.50-i-0.37N2-0.15/(N2-1-0.10) (5) 


and max(12-l-log(0/H)N2) = 8.86. For comparison. Pettini| 
& Pagel| ([2004jl linear calibration has parameters a = 8.90, 
b = 0.57, c = d = 0. The comparison of new N2 metallicities 
with 03N2 is shown in Figure® (right). Systematic offsets 
are removed. In the remainder of the paper we will only use 
recalibrated N2 metallicities. 

We return to the M*-Z relation, but now show the plot based 
on recalibrated N2 (Figure [T] right panel). SDSS loci of N2 
and 03N2 metallicities now agree better. As expected, there 
is less change in trends of z ^ 2 samples, except that the total 
change in N2 metallicity from low-mass end to high-mass end 
is now greater, i.e., the M*-Z trend is steeper. 

We now focus on high-redshift M*-Z trends in comparison 
with the local. The most striking difference is that the average 
N2 line ratios of the most massive galaxies (logM* > 11.0) 
are similar for local and high redshift galaxies (albeit there 
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Figure 3. |Salim et al.H2014) non-parametric analysis technique for studying the M,-Z-SFR relation. Metallicities of local and high-redshift galaxies are shown 
as a function of the relative SSFR (SSFR offset from the local star-for ming sequence), in four 0.5 dex-wide mass bins. Metallicities are derived using two 
methods: 03N2 metallicities (lower panels) use |Pettini & Pag^j2004) calibration, while N2 metallicities (upper panels) are based on our recalibration that 
forces agreement with 03N 2 for local galaxies (Figure]^. Kecalibration also caps N2 metallicities when N2 saturates. Black dots represent z ~ 2.3 galaxies from 
KBSS JSteidel et al.|2014), in compaiison with SDSS galaxies (grey points). High-redshift data show a statistically significant anti-con'elation at logM* = 10.0 
and 10.5 (panels B and (J), as indicated by the red line (linear weighted fit). Formal slopes and their errors are indicated, but whether they are actually statistically 
significant given the sample size requires additional consideration (Section 4.2). Visually the correlations appear less robust because the points with larger 
metallicity en'ors contribute to the scatter. Dashed lines show unweighted linear fits. Green lines are binned averages for SDSS galaxies (0.2 dex wide bins). For 
both indicators, SDSS and KBSS trends agree in the lowest mass bin, but become offset in subsequent bins, with 03N2 offset being somewhat larger than N2 
offset. 


is a large scatter), implying little metallicity evolution at the 
highest end, while the local and z ^ 2 M»-Z relations are quite 
offset in the case of 03N2, suggesting a strong evolution 
0.25 dex). In Sections 4.4 and 5.3 we will discuss possible 
causes of this relative discrepancy. 

Finally, from Figure [T] we also see that, regardless of the 
indicator, the metallicities of z ~ 2 galaxies do not saturate 
at high mass es as they do locally, as pointed out in |SteideI| 
|et al. ( |2014| l. KBSS and MOSDEF M*-Z relations gener- 
allyagreeror a given indicator, but with some differences 
in details, which will be discussed in Section 5.2. The av¬ 
erage trends of KBSS and MOSDEF samples are based on 
galaxies that are individually detected in requisite lines (Sec¬ 
tion 2). For MOSDEF, we do not show the average M»-Z 
trends below logM* ^ 10, where many individual galaxies 
are not dete cted in [NII]6584 line ( jKriek et al.|2014[|Sanders| 
|et al.|2015j l. Otherwise, KBSS and MOSDEF M*-Z relations 
based on indivi dual galaxies agree with the relations bas ed on 
stacked spectra (|Steidel et al.|2014| Sanders et al.|2015 1. 


4.2. Specific SFR as a second parameter of z ~ 2 
mass-metallicity relation 

M*-Z-SFR relation may be present at z ~ 2 even if it 
does not follow the local M*-Z-SFR relation (i.e., is not red- 
shift invariant), and we therefore hrst focus on establishing 
whether (S)SFR is a second parameter in KBSS sample. Fig- 
urel^applies the analysis framew ork for investigating the M*- 
Z-;^R relation, as introduced in |Salim et al.| ( |2OT4l l. N2 and 
03N2 metallicities are examined against the relative SSFR in 
four 0.5 dex wide mass bins, for both the SDSS and KBSS 
samples. SSFRs for both sampels are re lative to typical local 
values at that mass. [^03N2 is based on |Pettini & Pagel|2004 


^ We exclude from analysis Q2343-BX231, which, with calculated SFR = 


calibration, while N2 comes from our recalibration (Equation 

01 . 

To test for SFR dependence we perform weighted and un¬ 
weighted least square fits to KBSS points in each mass bin 
(full and dashed red lines in Figure]^. Weighted fits produce 
statistically signihcant non-zero slopes k (Z-SSFR correla¬ 
tions) in all mass bins for both metallicity indicators, except 
for logM* = 11.0 bin for N2 (k = -0.21 ±0.09). The val¬ 
ues of slopes are given in Figure In general, the slopes 
of KBSS galaxies are comparable to those of SDSS galax¬ 
ies at the same mass and SSFR. More specihcally, like SDSS 
galaxies, KBSS sample shows an anti-correlation between Z 
and SSFR in logM* = 10.0 and 10.5 bins (k = -0.14 ±0.02 
and -0.13 ±0.02 for N2 and -0.17 ±0.01 and -0.13 ±0.02 
for 03N2), for the two bins respectively. 

The formal slope errors do not take into account the uncer¬ 
tainties arising from small sample size. Therefore, to assess 
the statistical signihcance of anti-correlations we perform two 
additional tests: (a) we refit 100,000 bootstrapped samples 
and (b) we perturb the measurements by randomly drawing 
from a gaussian with cr equal to the reported metallicity error 
and by 0.15 dex for SSFR, and reht the trends 100,000 times. 
In both tests we perform weighted hts. At logM* = 10.0, the 
bootstrap test yields an anti-correlation in 94% of cases for N2 
(99% for 03N2), while perturbed fits have a negative correla¬ 
tion in 99% of cases (100% for 03N2). Results are similar for 
logM* = 10.5 bin: anti-correlation is present in 90% (96% for 
03N2) of bootstraps and 91% (93% for 03N2) of perturbed 
hts. In the highest mass bin, the trend of KBSS galaxies is 
not statistically signihcant (there is a similar number of cor¬ 
related and anti-correlated bootstraps). Currently the sample 

500 M0yr“* is an outlier in z ~ 2.3 SSFR—M*relation, lying 1.4 dex above 
it, and is also an outlier in the Z-SSFR relation. 
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Figure 4. Comparison of 03N2 and N2 metallicities of high-redshift samples (KBSS, left; MOSDEF, right) and local “analogs” to KBSS (middle). We define 
an analog to be a SPSS galaxy with the same SSFR and stellar mass as a KBSS galaxy. N2 metallicity is based on our recalibration which forces agreement with 
IPettini & PagelH2004I 03N2 metallicities (Figure|2) for typical local galaxies (grey points). Local analogs to KBSS sample have a similar offset as KBSS and 
MUSUEF samples, validating the M,-Z-SFR analysis approach th at is based on relative c omparison of metallicities of high-redshift and local galaxies of the 
same mass and SSFR. Possible AGN (N2> —0.5 and lying above the|Kauffmann et al.|2003|line) are shown as cyan symbols and were excluded from calculating 
the average offset in N2 metallicity (A), shown in each panel. The irow m the left panel shows the direction and a possible magnitude of AGN contamination 
vector affecting the cyan points. 


size in this bin is too small to establish if the metallicities of 
z^2 galaxies, like the local galaxies of similar mass, lack the 
dependence on SSFR. In the lowest mass bin KBSS galax¬ 
ies do not show an anti-correlation with SSFR, but neither do 
SDSS galaxies with such masses and SSFRs. While the over¬ 
all SSFR dependence in SDSS at logM* = 9.5 is quite strong, 
it appears to reach a low-metallicity plateau for SSFRs 1 dex 
or more above the main sequence. KBSS galaxies, which have 
such high SSFRs, also appear to have nearly contant metallic¬ 
ities. 

We note that the visual impression of the signihcance of 
some of the correlations differs from the results of the formal 
analysis that, unlike by-eye estimate, takes metallicity errors 
into account. Indeed, the significance of slopes having a non¬ 
zero value when fitting is performed without weighting by 
metallicity error (dashed lines in Figure is low, except for 
03N2 galaxies in logM* = 10.5 bin (k = -0.16 ±0.06). For 
the same reason (ignorance of weights), the visual impression 
of the vertical position of the best-fit linear trend does not al¬ 
ways agree with the its position based on the weighted fit (e.g., 
logM* = 10.0 bin for N2). 

Overall, we conclude that SSFR does appear to be a sec¬ 
ond parameter in z ~ 2 M*-Z relation, at least in the range of 
masses where such dependence is clearly present in the local 
relation. 

4.3. Redshift invariance of the mass-metallicity-SFR 
relation: standard assumptions 

We now turn our focus to the question of the invariance 
of M*-Z-SFR relation. In this section we will approach the 
analysis with standard assumptions regarding the local com¬ 
parison sample and the interpretation of high-redshift line ra¬ 
tios. In subsequent sections we will investigate the effects of 
possible AGN contamination of high-redshift line ratios, and 
of high-redshift selection effects due to high [OIII] sensitivity 
threshold. 

If the M»-Z-SFR relation is redshift invariant, then at any 
given mass and SSFR, the average metallicities of SDSS 
galaxies should be statistically consistent with those from 
KBSS. Therefore, Figure displays N2 and 03N2 line ra¬ 
tios as metallicities. However, significant discrepancies be¬ 
tween N2 and 03N2 metallicities of high-redshift galaxies 
were found when local calibrations were used to derive them 


([Newman et al.|2014HZahid et al.|2014bHSteidel et al.|2014[ 

Sanders et al.|2015|l, suggesting “evolution” in one or both in- 
dicators. These offsets can be seen in Figure [^for KBSS (left) 
and MOSDEF (right) samples. It may thus appear that possi¬ 
ble “evolution” of metallicity calibrations would preclude the 
test of redshift invariance of an M»-Z-SFR relation. How¬ 
ever, Figure (middle panel) shows that, remarkably, simi¬ 
lar offset between N2 and 03N2-inferred metallicities is also 
present in local high-SSFR galaxies (we specifically show 
what we call the local “analogs” of the KBSS sample, see 
Section 3). We interpret this to mean that the critical limita¬ 
tion for the application of local calibrations at high redshift 
is not so much that they are local, but rather that they are 
based on typical local galaxies (or HII regions in such galax¬ 
ies). Therefore, the relative comparison of line ratios of high- 
redshift galaxies and local galaxies of similar SSFR, which 
is at the essence of establishing the invariance of any poten¬ 
tial M*-Z-SFR relation, can be carried out regardless of the 
uncertainties involving the conversion of line ratios into abso¬ 
lute metallicities, i.e., without having to decide if it is N2 or 
03N2 (or both) indicator for which the “local” calibration is 
off at high redshift {and for local galaxies with high relative 
SSFRs). With this important conclusion in hand, we proceed 
with the analysis. 

Figure compares metallicities of SDSS (grey dots with 
binned averages shown with green line) and KBSS galaxies 
(black points with linear fits shown by red line), using N2 
(upper panels) and 03N2 (lower panels) metallicity indica¬ 
tors. We notice that SDSS galaxies extend into the range oc¬ 
cupied by KBSS, i.e., in each bin there do exist local galaxies 
with SSFRs as high as those in z ~ 2.3 sample. While such 
galaxies are extremely rare today, they are present in the large 
volume probed by SDSS. Hence, direct comparison of local 
and z ~ 2 populations is possible, and parameterizations of 
the local Mf,-Z-SFR relation do not need to be extrapolated 
into the range of SSFRs occupied by high-redshift galaxies in 
order to test for possible evolution. 

Examination of Figure shows that the average trends of 
both N2 and 03N2 line ratios for SDSS and KBSS gener¬ 
ally agree at lower masses. However, an offset appears in 
logM* = 10.0 bin, in the sense that KBSS metallicities are 
lower, especially for 03N2. At logM* = 10.5 the offset grows 
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Figure 5. Metallicity difference of KBSS galaxies and their local “analogs” 
in SDSS using N2 (upper panel) and 03N2 (lower panel) indicators. Analog 
is a galaxy with matching SSFR and stellar mass. The offset is mass depen¬ 
dent. The vertical eiTor bar incl udes only the uncertainty of KBSS metallicity, 
as listed in|Steidel et al.|(2014), so it is a lower limit of the full uncertainty of 
the metallicity difference. Colored line show binned averages. Dashed lines 
show the interval centered on logM* = 10.5 used to calculate the average 
offset at high mass (A 10 . 5 ). 

to ^ 0.25 in 03N2, but is somewhat smaller for N2, where it is 
compensated by scatter towards higher metallicities. The off¬ 
set stays large for 03N2 at logM* = 11.0, but is again some¬ 
what smaller for N2. Interestingly, even in bins in which the 
offsets between average trends are large, there typically do ex¬ 
ist SDSS galaxies with such low metallicities as KBSS galax¬ 
ies. 

The exact degree of offsets between the SDSS and KBSS 
trends depends on how one chooses to treat the data. If metal¬ 
licity weights are ignored, as in fitting by-eye, the trend for, 
e.g., N2 in logM* = 10.0 bin becomes more offset from the 
SDSS (dashed line). Thus we also perform an additional, 
more direct test of invariance. Instead of comparing the trends 
in various mass bins, we now study the difference in metal¬ 
licity between KBSS galaxies and their local “analogs” from 
SDSS (i.e., the galaxy with matching SSFR and M* , see Sec¬ 
tion 3, and lying below the Kauffmann et al.|2003| AGN de¬ 
marcation line). We show the differences in metallicities for 
each KBSS-analog pair against its stellar mass (Figure |^. 
Redshift-invariant relation requires a lack of systematic OTt- 
sets. We see that this is the case below logM* ^ 9.7, but 
above it the systematic difference appears, confirming the re¬ 
sults from Z-SSFR analysis. In a 0.5 dex interval centered at 
logM* = 10.5 the average difference is -0.23 ± 0.04 dex for 
N2 and -0.28 ± 0.03 dex for 03N2 (errors of offsets were ob¬ 
tained from bootstrapping resampling; they are not standard 
deviations). 

4.4. Redshift invariance: AGN contamination 

From the analysis presented so far we would conclude that 
the M*-Z-SFR relation is altogether not invariant, because 
of the systematic differences in metallicities at larger masses. 


Figure 6. Metallicity difference of KBSS galaxies and their local “analogs” 
in SDSS. Same as Figure]^ except that we now apply a correction to galaxies 
whose line measurements may be affected by an AGN contribution (cyan dots 
in Figurel^ left). The overall effect of the correction is relatively subtle, but 
it does produce more consistent offsets for N2 and 03N2 indicators, for the 
entire sample (A), and at higher mass (Aio.s). 


The differences tend to be larger using the 03N2 than N2 in¬ 
dicator. This discrepancy could be the result of unrecognized 
AGN contribution to emission lines. Namely, AGN contribu¬ 
tion moves a galaxy along the AG N mixing sequence in the 
BPT diagram Baldwin et al. (|1981|l (i.e., 03 vs. N2 diagram). 


increasing both N2 and 03. The increase in 03 is larger than 
in N2, so net 03N2 ratio increases as well. Larger N2 gives 
higher apparent metallicity, while larger 03N2 gives lower 
apparent metallicity, giving rise to a discrepancy. 

In Figure|7](lower left) we show the BPT diagram for KBSS 
sample. Some galaxies, especially with larger N2 values 
(N2> -0.5) are found away from the locus of the rest of the 
galaxies, along what may be an AGN spur. These galaxies 
tend to be more massive. These possible AGN-contaminated 
galaxies are shown as cyan dots in Figure|^(left), in which we 
directly compared derived 03N2 and N2 metallicities. These 
candidate AGN indeed lie even further from the 1:1 diagonal 
th an the rest of the KBSS sample. Similar effect was reported 
in|Newman et al.|(|2014|l. Moving a galaxy up the AGN mix¬ 
ing sequence increases 03 and N2 such that A03 « 5AN2. 
For A03 = 0.4 dex, which appears as a reasonable value 
from inspecting Figure the resulting 03N2 “metallicity” 
decreases by 0.10 dex, ^ile the N2 “metallicity” increases 
by 0.12 dex, i.e., galaxy is shifted almost perpendicular to the 
diagonal in Figureffl(arrow in the left panel). Similar shift is 
seen in several MO^EF galaxies (right panel). 

Is it realistic to have unrecognized AGNs in the KBSS (and 
MOSDEF) samples? KBSS excluded high ionization species 
emission-line galaxies as being contaminated by AGN. This 
method is only sensitive to more energetic AGNs (ionization 
potentials for C IV and N V are 48 and 77 eV, respectively, 
compared to 15 eV and 35 eV for [Nil] and [OIII]). A similar 
method applied to the SDSS would not identify as AGN many 
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Figure 7. Comparison of z ~ 2.3 KBSS (left), their local SDSS analogs (middle) and z ~ 2.3 MOSDEF samples (right), in terms of SSFR-M* plane (upper 
row), and the BPT diagrams (lower row). Green lines in the upper row represent a linear fit to the SDSS SF sequence (Equation]^ greyscale). High-redshift 
samples, especially KBSS, have steep slopes of SSFR vs. mass compared to local galaxies, or even compared to other surveys at that redshift (3D-F[ST,|Whitaker| 
|et al.|2014I. Fiigh-redshift surveys cannot detect emission lines of potential galaxies in the lower part of the BPT diagram (red and blue contours show SDSS 
galaxies with [OIII] and Ha luminosities above KBSS and MOSDEF detection limits). The lower part of the BPT diagram is where locally the more massive 
galaxies are found (yellow dots in left and middle lower panel have logM* > 10.25). Dashed vertical line in lower left panel shows separates on the right possible 
AGN in KBSS. Rectangle in the lower left panel displays the selection of SDSS galaxies that mimics th e region detected and occupied by KBSS sample, which 
admits some galaxies that lie above the local AGN demarcation line (blue curve,[Kauffmann et al.|2Q0^. 


of the galaxies on the AGN branch of the BPT diagram. Line 
ratio diagnostics may be the only practical method to recog¬ 
nize weaker AGN even locally. 

How does the possible AGN contamination bear on test¬ 
ing the M*-Z-SFR redshift invariance? In Figure|^we again 
show the difference between the metallicities of KBSS galax¬ 
ies and their local analogs as a function of mass, but we 
now correct the metallicities of AGN candidates by -0.10 
dex for N2 metallicity (13 galaxies) and -1-0.12 dex for 03N2 
metallicity (11 galaxies). Systematic offsets remain, espe¬ 
cially at high masses, but they are now nearly identical for 
N2 and 03N2 metallicities (-0.13 dex overall and -0.26 dex 
at logM* = 10.5). 

4.5. Redshift invariance: [OIII] sensitivity threshold 

In this section we explore w hether the sensitivi ty of high- 
redshift spectroscopic surveys Puneau et al.|2014| affects the 
analysis of the evolution of the M*-Z-SFR relation. 

To illustrate the potential issue with sensitivity, we again ex¬ 
amine the BPT diagrams in Figurel^(KBSS, lower left panel; 
MOSDEF, lower right panel). Botn z ^2.3 samples are con¬ 
fined to the upper part of the diagram, which is locally domi¬ 
nated by lower metallicity galaxies. What is interesting is that 
the regions of the BPT diagram to which the current high- 
redshift spectroscopic surveys are sensitive to (shown with 
red and blue contours in KBSS and MOSDEF BPT diagrams, 
respectively), largely coincide with the locus of individual de¬ 
tections. (The sensitivity contours were obtained by selecting 
SDSS galaxies having [OIII]5007 and Ha luminosities above 


KBSS/MOSDEF survey limits (Section 2), following |Juneau| 
et al.|2014 ) If galaxies at z ~ 2 existed in the lower portion of 


the BPT diagram, current surveys would not include them in 
their samples, because they require [OIII]5007 line detection. 

We can also use local “analogs” to illustrate the nature 
of this potential issue. Like KBSS and MOSDEF, the local 
“analogs’ll to KBSS galaxies are offset upwards in the BPT 
diagram (Figure 1^, lower middle panel), but unlike the z ~ 2 
samples, the locm “analogs” to the KBSS sample populate 
the full extent of the star-forming branch in the BPT diagram, 
with the majority of the more massive “analogs” having low 
03 (yellow points indicate logM* > 10.25). The difference in 
the extent of the locus of the “analogs” on the BPT diagram 
with respect to z ~ 2 samples is essentially another way of 
saying that the M*-Z-SFR relation shows no offset at lower 
masses (upper part of the BPT diagram), but becomes increas¬ 
ingly discrepant at higher masses, because the locally higher- 
mass galaxies populate the lower regions of the BPT diagram, 
while they are (intrinsically, or perhaps because of sensitivity 
limits) found only in the upper parts at high redshift. 

To understand the possible consequences of the insensitiv¬ 
ity of high-redshift surveys to galaxies with low 03 ratios 
(high metallicities), we repeat the analysis, but now select 
SDSS galaxies only from the region of the BPT diagram ac¬ 
cessible to, and occupied by KBSS galaxies (0.1 < 03 ^ 0.9, 
and -1.7 < N2 < -0.3, shown with the rectangle in Fig- 


In this figure we allow the “analogs” to be selected regardless of the 
position i n the BPT diagram, i.e., including from among the galaxies lying 
above the |Kauffmann et ar]j2QQ3| AGN demarcation line. 
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Figure 8. Dependence of N2 and 03N2 metallicity on the offset from the local star-forming sequence for z ~ 2.3 galaxies (black dots) and SDSS galaxies (grey 
points) selected from a similar region of the BPT diagram as occupied by high-redshift samples (rectangle in Figure!^ lower left). Selection is modified from 
the standard one in order to allow for the possibility that high-redshift surveys miss galaxies in regions of the BPT oTagram that they are not sensitive to, and 
to allow for the possibility of AGN contamination. The agreement between line ratio (metallicity) trends is improved, i.e., in this scenario the redshift invariant 
M*-Z-SFR relation would be viable. However, other evidence disfavors the existence of a z ~ 2 high-metallicity population. See also Figure[^caption, except 
that binned averages for SDSS are now shown as orange lines. 



Figure 9. Metallicity difference of KBSS galaxies and their local “analogs” 
in SDSS. Analogs are now selected from among the SDSS galaxies that oc¬ 
cupy a similar region of the BPT as high-redshift samples, in order to account 
for the possibility that high-redshift surveys miss galaxies in regions of the 
BPT diagram that they are not sensitive to, and to allow for AGN contami¬ 
nation. The offsets are reduced compared to standard selection/assumptions 
(Figure]^, but are not entirely eliminated. 

ure [7] lower left). This selection box will also admit some 
SDSS galaxies that are found above the AGN demarcation 
line, therefore implicitly “correcting” for the possibility that 
z ~ 2 samples are also affected by AGN contribution. We 
therefore do not apply any explicit AGN correction as we did 


in Section 4.4. 

The new set of Z-SSFR plots is presented in Figure!^ The 
offset between SDSS and KBSS trends has generally rMuced 
compared to FigureThe large offset in the logM* = 11 bin 
is eliminated, and is greatly reduced in the logM* = 10.0 and 
10.5 bins. The improved agreement arises primarily because 
the local galaxies with high metallicity, which the z ^ 2 sur¬ 
veys would not be sensitive to, are now excluded from the 
comparison. These galaxies pushed the average metallicity 
trend of SDSS galaxies upwards. 

We also repeat the analysis in which we directly compare 
the metallicities of KBSS galaxies with their local “analogs”, 
except that we select the analogs from the rectangle in Figure 
It] lower left. The average trends of metallicity difference in 
Figure 1^ confirm that the discrepancies at higher masses are 
significantly reduced, although they are not completely elimi¬ 
nated. At logM* ~ 10.5 the offset is now -0.05 ± 0.04 for N2 
and -0.08 ± 0.03 for 03N2. N2 and 03N2 offsets are compa¬ 
rable. The overall offset (at all masses) is now -0.06 ± 0.02 
with N2, and nearly identical with 03N2. 

The above analysis demonstrates that our ability to con¬ 
strain the metallicity evolution at higher masses hinges on an 
implicit assumption used in all studies so far that there do 
not exist significant populations of high-redshift galaxies in 
regions of the BPT diagram that are not accessible to obser¬ 
vations (low 03 values). May such high-redshift populations 
exist? The fact that the ISM conditions at z ~ 2 appear to be 
so mewhat more extr eme than they are in typical local galax¬ 
ies dCoil et al.|2015| l does not preclude populating the lower 
(high-metallicity) portion of the BPT diagram. This is evident 
from the position of local analogs in Figure (middle lower 
panel), which follow the “more extreme” location of KBSS 
galaxies in the upper part of the BPT diagram, but neverthe¬ 
less also populate the lower portions. It is obviously important 
to understand if observational threshold actually produce any 
such biases. Because we do not have access to original KBSS 
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and MOSDEF data, we base this discussion on the inferences 
fr om published results. 

i Shapley et al.| (|20151l tested whether non-detections create 
)ias in the BET distribution for the MOSDEE sample. They 
produced average-stacked spectra, in mass bins, that included 
both the detections and the non-detections. The stacked spec¬ 
tra had line ratios that placed them in the region of BPT values 
shared by the detections (03> 0.2, their Fig. 2, left). Because 
the line ratios are distributed log-normally, median stacking of 
spectra would be preferred to average stacking, however, the 
results of such exercise produce similar results (A. Shapley, 
priv. comm.) A more stringent test would involve stacking 
only the non-detections. However, even the stacked spectra 
of non-detections are unlikely to fall in the lower part of the 
BPT diagram because the majority of non-detections in both 
KBSS and MOSDEF are galaxies with weak [Nil], and not 
with weak [OIII]. Furthermore, the detection fraction tends 
to be high (> 80%) at higher masses. Therefore, it appears 
that high-redshift samples, modulo some other sample selec¬ 
tion bias (Section 5.1), are not biased by high [OIII] detec¬ 
tion thresholds, i.e., they are probably not missing a signif¬ 
icant population of high-metallicity galaxies. Consequently, 
the fact that the position of z 2 samples in the BPT diagram 
matches the regions that these surveys are sensitive to (Figure 
1^ appears to be a mere coincidence, and a gap between local 
and high-redshift M*-Z-SFR relations cannot be attributed to 
sensitivity issues. Nevertheless, we suggest that future studies 
should consider this issue. 


5. OPEN QUESTIONS 

5.1. Are KBSS and MOSDEF spectroscopic samples 
representative of star-forming galaxies at z ^ 2? 

In order to robustly assess the evolution of chemical en¬ 
richment from z ^ 2 to the present day, one requires both 
the high-redshift and the local samples to represent typical 
star-forming galaxies at that redshift. It is currently uncer¬ 
tain whether this is the case for high-redshift spectroscopic 
samples. The distribution of KBSS and MOSDEF data in 
SSFR-M* plane is different from what is expected. The slope 
of the star-forming sequence of KBSS sample is very steep 
(/3 = -0.84), signihcantly steeper than /3 = -0.35 slope of the 
local galaxies (Equation 2; green line in the upper panels of 
Figure]^. This disagrees with most results that suggest that 
the slope at z 2 is similar or s hallower than the local slope 
( [Speagle et al.|20f4 i. Recently, [Whitaker et al.| ( [2014) l have 
determined the star-forming sequence at z 2.3 based on an 
extensive dataset from 3D-HST survey (purple line in Fig¬ 
ure]^ upper panels). As expected, the slope of the 3D-HST 
main sequence is somewhat shallower than the local one (we 
get P = -0.19), in sharp contrast to KBSS. At logM* < 9.7 
the 3D-HST sequence is lower than the KBSS sequence (but 
not MOSDEF), probably due to the KBSS’s UV selection, as 
mentioned in jSteidel et al.| ( |70141 l. However, the situation re¬ 
verses at logM, > 9.7, where both KBSS and MOSDEF have 
lower SSFRs than 3D-HST, up to 0.4 dex at logM* « 11. This 
high-mass offset has not been discussed in the literature so 
far. Altogether, it appears that the z 2 spectroscopic sam¬ 
ples target somewhat more intense star-formers at lower mass 
and more quiescent galaxies at higher mass. This potential 
bias does not directly affect for our M*-Z-SFR relation anal¬ 
ysis because we are comparing SDSS and KBSS galaxies at 
the same specific SFR. Nevertheless, it is important to fully 
understand if spectroscopic samples are representative of un¬ 


derlying star-forming population at z ~ 2, especially at high 
mass, where we hnd that the M*-Z-SFR relations have an 
offset. 

5.2. Are there systematic differences between KBSS and 
MOSDEF line ratios? 

The M^-Z relations produced using KBSS and MOSDEF 
data on individual detections agree very well (Figure [T]), ex¬ 
cept above logM* = 10.2 for 03N2, where MOSDEF metal- 
licities start to diverge from KBSS ones, the latter being lower. 

The discrepancy between MOSDEF and KBS S measure¬ 
ments has been discussed in Shapley et al. P015 i in terms of 
the position of the galaxies in the upper portion of the star¬ 
forming branch in the BPT diagram. Here we perform similar 
analysis, but quantify the differences between the samples by 
hnding the value 03med, such that the curve: 


03: 


0.61 


N2-1-0.08 


-|-03nied 


(6) 


divides the sample in half. This function has only one free 
parameter: 03ined- The values of the other two para meters 
are hxed to the values used by Kewley et al. ([2013 1 to de 


scribe the star-forming locus at z = 0. jSteidel et al.j ([2014 1 
have allowed all three parameters to be free, but we niid that 
the single-parameter form with other parameters hxed to lo¬ 
cal values actually better describes the high-redshift samples 
after possible AGN-contaminated galaxies are taken out. Fur¬ 
thermore, a single-parameter expression, i.e., a shift in the 
03 direction, ma kes it easier to comp are different samples. 
For local galaxies [Kewley et al. P013|l hnd 03med=l l- For 


high redshift samples and their local analogs we focus only on 
the upper left part of the BPT diagram (N2< -0.5, 03> 0.0), 
thus excluding possible AGN contamination (see also Figure 
0 and outliers. For KBSS we obtain 03ined=l-32, i.e., its star¬ 
forming branch lies some 0.2 dex higher in the BPT diagram 
than the local one. For MOSDEF, we get O3[ned=l -20, which, 
while higher than the local value, is signihcantly lower than 
the KBSS value. Local analogs to KBSS yield 03med=l-28, 
which is closer to KBSS than MOSDEF, and also much higher 
than the typical local galaxies. That the local galaxies with 
high relative SSFRs were offset in the BPT diagram was pre- 


viously found in 

Brinchmann et al. (|2008|l. 

Shapley et al. ( 

2015 1 ascribe the difference between the po- 


sition of KBSS and MOSDEF samples in the BPT diagram 
sample selection: UV selection of KBSS sample as opposed 
to the optical selection of MOSDEF. However, it is not clear 
whether selection is the main culprit. Inspecting the SSFR- 
M* distributions of the two samples (Figure]^ upper left and 
right), they appear similar, except at logM* < 9.5 where 
KBSS sample has higher SSFRs and extends to lower masses. 
Even when we restrict the KBSS sample to mass range that 
both KBSS and MOSDEF appear to cover in similar fashion 
(9.5 < logM* < 10.5), we still obtain 03med=l-29, a signifi¬ 
cantly higher value than for M0SDEF[^ Local analogs in the 
same mass range yield 03med=l-27, i.e., very close to KBSS 
value. Future studies should be able to pin down the source of 
the discrepancy. The analysis presented here suggests that the 
cause may not be solely the differences in sample selection. 

We also wish to point out that while 03[ned are similar for 
KBSS and their local analogs, the former have a larger scatter 
(presumably because of larger measurement errors), which is 

** We cannot perform equivalent, mass-range limited determination for 
MOSDEF because we do not have matched tables of line ratios and masses. 
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Figure 10. Evolution of the mass-metallicity relation from z ~ 2 to today. 
Concordant high-redshift trend (violet band) is a qualitative illustration 
of how to reconcile the M*-Z relations based on 03N2 (solid lines) and N2 
(dashed lines) metallicity indicators, from KBSS (red) and MOSDEF (blue) 
surveys. To do so we consider various systematics discussed in Section 5.3 
and labeled in the figure. Generally, at lower masses the concordant trend is 
closer to 03N2 because N2 will be more affected by the change in the N/0 
ratio, and at higher masses it is between N2 and 03N2 trends because of the 
possible AGN contamination that boosts N2 but suppresses 03N2. 


why they more often cross the local AGN demarcation line 
than their analogs, which are found close to the line but do 
not cross it. Indeed, if we increase 03 values of analogs 
by the difference in 03med between the analogs and KBSS 
(0.046 dex), only 6 out of 71 galaxies fro m the upper left of 
the BPT diagram would actually cross the |Kauffmann et al.| 
(20031 line. This suggests that z ~ 2 galaxies without AGN 


contribution intrinsically mostly lie below the Kauffmann et 


al. (|2003 1 line. This conclusion would be even stronger for 


MOSDEF, whose 03med is lower than that of KBSS analogs. 
We conclude that ha d the z ~ 2 line ratios had the same preci¬ 
sion as in SDSS, the [Kauffmann et ^^003) 1 AGN classifica¬ 
tion line would not require a shift greater than 0.05 dex in the 
vertical direction. The appearances are significantly differ¬ 
ent because of the larger measurement error that contributes 
to the scatter and because some of the galaxies that are cur¬ 
rently not considered to be AGN probably are AGN. Allowing 
for larger observational errors, but still removing the potential 
AGN may be achieved using a simple cut in 03 and N2: 


N2>-0.5andO3>0.1, (7) 

which resonates with the conclusions of [Coil et al.| ( |2015T ) that 
were based on the MOSDEF sample. 

5.3. Towards the concordant evolution of mass-metallicity 

relation 

In this section we will attempt to arrive at a consistent M*- 
Z relation at z 2 and discuss its evolution. This requires 
the knowledge of absolute metallicities. As discussed in Sec¬ 
tion 4.3, the analysis of M»-Z-SFR invariance, being relative, 
could be carried out even without the knowledge of absolute 
calibrations that would be needed to convert the line ratios of 
high-redshift galaxies (and their local analogs) into metallici¬ 
ties. This is not the case if we are interested in constructing the 


M*-Z relations. The remaining offset between N2 and 03N2 
“metallicities” eve n after we have correc ted for a mismatch 
that was present in Pettini & Pagel ( |2004| l relations, and have 
corrected for possible AGN contribution (Figure |^, means 
that either or both of the local calibrations does not hold for 
ISM conditions in z ~ 2 galaxies and their local analogs. 

Whether N2 or 03N2 is more affected by the ch anges in 
the ISM conditions is currently a matter of debate. [Masters 
et al. ( 2014[ ) and jShapley et al. ( 2015[ l argue for 03N2 be¬ 
ing less affected because z 2 galaxies appear to follow lo¬ 
cal galaxies in 03 vs. S2 (= log([SII]6573/Ha)) and 032 vs. 
R23 diagrams. Such diagnostics may not be particularly rel¬ 
evant in the context of 03N2 vs. N2 debate because using 
03, and especially 032, results in comparing high-redshift 
galaxies with local galaxies having high ionization parame¬ 
ters, which, as we have shown, also exhibit the 03N2 vs. 
N2 discrepancy. [Steidel et al.[ ([2014[) consider 03N2 to be 
less biased than N2 based on a very small number of direct 
metallicities of the local “green pea” galaxies, while[Liu et al. 
( [2008 1 reach a similar conclusion based on direct metallici- 
ties of stac ked SDSS galaxies w ith high 03 values. On the 
other hand, [Zahid et al.[ ( [2014b| l consider N2 to be more reli¬ 
able based on the greater sensitivity of 03N 2 to changes in the 
ionization parameter, which many studies (Brinchmann et al. 


2008[ Nakajima & Ouchi|2014| [Kewley et al.[[2013p assume 

to be the principal difference between IS M conditions of lo- 
cal and high-redshift galaxies. In contrast, Steidel et al.|2014 
show, using photoionization models, that the ionization pa- 
rameter cannot be the only factor driving the difference. The 
greater effect of the ionization parameter on 03N2 may be 
countered by higher N/O ratio at z 2 compared to typical 
lo cal galaxies of t he sam e metallicity, as tentatively measured 
by Masters et al. ( 2014[ l, thus making 03N2 more reliable for 
high-redshift stuies than N2 in the end. 

For the purposes of the remaining discussion we will as¬ 
sume that 03N2 provides more robust metallicities than N2 
and that the local 03N2 calibration can be applied at z 2. 
Figure 0 then implies that to correct our recalibrated N2 
metallicities for high-redshift/local analog galaxies one would 
need to subtract ~ 0.08 dex. 

Based on all of the inferences made so far we present a 
possible evolution of the M*-Z relation from z 2 to today in 
Figure!^ This figure shows KBSS and MOSDEF metallicity 
trends based on both N2 (dashed lines) and 03N2 (solid lines) 
and then selects a trend that better fits available evidence (vi¬ 
olet band). At lower masses (logM* < 10.3) the final trend is 
closer to 03N2 than to N2 because 03N2 indicator is likely 
more robust for samples with very high SSFRs. We allow the 
final trend to sit slightly above the 03N2 trends to account 
for z ~ 2.3 samples possibly having atypically high SSFR for 
that redshift (Section 5.1). As we progress towards higher 
masses, the possible AGN contamination affects both 03N2 
and N2, but in different directions. Applying a AO3=0.4 dex 
correction to AGN candidates (cyan dots in Figure results 
in N2 and 03N2 trends that are much closer to each other and 
roughly between the uncorrected values. 

Comparing the final z 2 trend now with the local M*- 
Z relation (green band) we see that the level of metallicity 
evolution appears to be mass-dependent above logM* = 10.3, 
because unlike the local M*-Z relation, the high-redshift one 
does not saturate (flatten). Note that the proper comparison 
requires the local M*-Z relation to be constructed in an un¬ 
biased way. In particular it is important not to introduce an 
[OIII]5007 selection, as it would preferentially remove the 
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highest metallicity galaxies and could therefore modify the 
character of the M^^-Z relation at higher masses. At higher 
masses the final z ^ 2 trend approaches that of the local galax¬ 
ies, but the gap is not fully closed. At masses below the local 
plateau, the z ~ 0 and z 2 M*-Z relations have a similar 
sl ope, as advocated in “universal metallicity relation” model 
of |Zahid et al.| ( [2014a| l. 

6 . DISCUSSION 

6.1. Prior work on z ~ 2 M^-Z-SFR relation 

[Steidel et al.| ( |2014| ) and [Sanders et al.| ( |2015| l did not de¬ 
tect a significant secondary dependence of the M*-Z relation 
on SFR in their analysis of KBSS and MOSDEF data, while 
our analysis of some of these same data does find a depen¬ 
dence on SFR. The likely reason why the SFR dependence 
was not found in previous z ~ 2 studies is because the anal¬ 
ysis methods used in those studies did not fully account for 
the fact that the M»-Z-SFR relation is intrinsically not very 
tight, even locally where measurement errors are smaller, and 
that the relation is driven by changes in SFR relative to the 
typical SFR at a fixed mass , rather tha n absolute SFR ( Salimj 
jet al.|[2OT4| ). [Steidel et al.| ( |2014l ) and [Sanders et al.| ( 2015[ l 
split the sample into high and low SFR and then look for off- 
sets between respective MZRs. This method will not show a 
SFR dependence because it selects by absolute SFR, and be¬ 
cause the information from a limited range of SFRs at a given 
mas s is co llapsed into two closely separated bins. [Wuyts et] 
[ak] ( |2014| ) split their sample into mass-dependent low/high 
Wr bins (also mass-dependent SFR quartiles, R. Sanders, 
priv. comm.), which is more similar to our methodology, but 
even the mass-dependent binning is apparently too crude to 
uncover a relatively weak SFR dependen ce, especially in rela¬ 
tively small samples. [Maier et ak ( [2014| l do not bin their data, 
so they tentatively detect the dependence on SFR in their sam¬ 
ple of 20 galaxies at z ^ 2.3. The key to being able to tease 
out the SFR dependence in relatively noisy high-redshift data 
is not to bin by SFR or any other SFR-related quantity, but to 
directly look at metallicity as a function of SSFR (or relative 
SSFR) and do so separately for galaxies of different stellar 
masses. 

As pointed out, th e “direct" approach used here and in 
T] ([2015J for testing redshift invariance of the 
ilation ha 


Sanders et al 
M.-Z-SFR re 


has an advantage over previously applied 
methods because it does not rely on extrapolations of param- 
eterizations which approximate the shape of the surface de¬ 
fined by local galaxies, which do not adequately capture the 
behavior of outlier populations such as galaxies with high SS- 
FRs, and can le ad to differing pred ictions for metallicities 
at high redshift ([Maier et ak||2014[l. Here we have shown 
that extrapolations of the local trends are not needed - di¬ 
rect comparison to local galaxies that occupy the same part 
of SSFR-M* space as z ~ 2.3 galaxies is possib le. The same 
conclusion was reache d in[Sanders et ak[ ( 2015[l, while some 
previous studies (e.g., Maier et ak|2014[ Steidel et ak|2014) 
considered the need to extrapolate the parametenzations of 
local trends to be essential. However, they ba sed that need 
on the ra nge of SDSS SFR values reported in [Mannucci et| 
[ak] ( [2010] ), who used fiber SFRs, which, in addition to being 
oistance- dependent, are on a verage 0.6 dex lower than the to¬ 
tal SFRs ( Salim et ak||2014| . Many studies have in addition 
adopted Mannucci et ak parametrization of the local trends 


studies, Sanders et ak (20151, like our study, avoided parame- 
terization, and used the line ratios from the actual SDSS data, 
binned by mass and total SFR, to perform a direct compari¬ 
son with their MOSDEF sample. They find that the z 2.3 
metallicities are ^ Q.ldex lower than SDSS metallicities of 
galaxies with similar mass and SFR, similar to our overall re¬ 
sult (Figure]^. We present a more nuanced picture where the 
discrepancy is preferentially present in more massive galax¬ 
ies. 

Our study also sheds new light on the discrepancy between 


N2 and 03N2 metallicities at z 2 ([Newman et ak 
Cullen et ak||2014[ [Zahid et ak|[2014b[ [Steidel et ak 


2014 


2014 


Sanders et ak|2015) . We find evidence for three independent 
causes. A smaller part (~ 0.03 dex) is due to a mismatch in 
local calibrations (Section 4.1). The larger part (0.06-0.09 
dex) is due to changed ISM conditions with respect to typical 
local galaxies. Figure ffldemonstrates that such different ISM 
conditions are not exclusive to high-redshift galaxies, but are 
common to their local analogs. Lastly, and preferentially af¬ 
fecting high-mass galaxies, an additional offset is due to AGN 
contamination of [Nil] and [OIII] lines, which increases the 
relative discrepancy in derived “metallicities” by additional 
^ 0.2 dex (Section 4.4). 

6.2. Implications 

Secondary dependence of the Mf,-Z relation on SFR at any 
redshift appears as a natural outcom e in recent models of 
galaxy evolu tion, both analytical ( e.g., [Lilly et ak|2013[[Dav^ 
[et ak [20T^ and numerical (e.g., [Dave et ak||2011| l It sig- 
nifies departure from equilibrium metallicity at a given mass 
due to the variations in the gas infall rate, which also mod¬ 
ulates star formation. That we should find evidence for SFR 
dependence at z 2 is therefore somewhat expected. How¬ 
ever, our findings show that discussions of SFR dependence 
must be more nuanced. The strength of the SFR dependence 
in the local universe obviously de pends on the mass, and is 
very weak or absent at high mass ([Ellison et ak||2008[ Man 


nucci et ak[2010[|Salim et ak|2014[|. This may be because the 

local massive galaxies (l ogM^, » 11), have by now reached a 
“saturation” metallicity ( [Zahid et ak|2014a) l, so are more dif¬ 
ficult to perturb chemically by gas infall, even though the in¬ 
fall would increase the SFR. Current samples at z 2 are too 
small to establish if high-redshift, high-mass galaxies behave 
in a similar way. They may not, considering that their metal¬ 
licities appear to be below the saturation level. The results we 
present here reveal that the situation may also be more com¬ 
plex at lower masses (logM, « 9.5). There, the overall local 
SFR dependence is very strong, but then appears to reach a 
low-metallicity plateau, such that the further increase in SSFR 
does not lead to further decrease in metallicity. High-redshift 
galaxies with logM* ~ 9.5, which have similar SSFRs as the 
local plateau galaxies also show no anti-correlation between 
metallicity and SSFR. Future theoretical studies should at¬ 
tempt to explain these detailed behaviors. 

For the M*-Z-SFR relation to be invariant with redshift re- 
quires addi t ional c onstraints. In the “gas regulator” model of 
Lilly et ak ( 2013|l, an invariant M*-Z-SFR relation emerges 


to test for M^.-Z-SFR invariance (e.g ., Cullen et ak 2014 
[Wuyts et ak||2014^ [Zahid et ak||2014bl l. In contrast to these 


only if the gas consumption timescale and m ass loading of 
wind outflows are const ant in time (see also [Forbes et ak| 
|2014[l. [Lilly et ak|([2013|l models with constant SF efficiency 
also predict that the metallicity evolution will decrease with 
mass, which, in the context of invariant M*-Z-SFR relation, 
is equivalent to having the SSFR dependence of M*-Z-SFR 
relation decrease with mass. Mass-dependent evolution of the 
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M^-Z relat ion is also predicted by the “universal metallicity 
relation” of|Zahid et al.|(|2014a|l, which relates metallicity and 
gas-richness. More detailed overv iew of recent theoretical ef¬ 
forts is given in |Salim et al.| ( |2014) l. Establishing whether M*- 
Z-SER relation is mvanant or not is therefore important to 
guide our understanding of the complex interplay between in¬ 
fall, outflows and SE. Again, our results challenge simple ex¬ 
planations. We find that local and z ~ 2 M*-Z-SER relations 
are consistent with each other at lower masses (logM* < 10), 
but then quickly reach a significant offset (~ 0.25). This re¬ 
sult suggests that the low-mass local “analogs”, rare galaxies 
that are found 0.7 dex or more above the main sequence, may 
have a similar mode of SF as the high-redshift galaxies of the 
same mass, and are similarly “unevolved”. High mass galax¬ 
ies, on the other hand, tend to be metal-rich today, even when 
their SFRs are as high as those of high-redshift galaxies, sug¬ 
gesting that they owe high SFRs to different processes from 
those that operate at high redshift. 


7. CONCLUSIONS 

We conclude that barring some selection effect, the M*-Z- 
SFR relation of z ~ 2 galaxies is consistent with the local one 
at lower masses, but not at higher masses, so it is overall not 
redshift invariant. We summarize the main conclusions of this 
study are as follows: 


1. Local (z ~ 0) galaxies with SSFR and M* typical of 
z ~ 2.3 spectroscopic samples exist in SDSS, which al¬ 
lows the redshift invariance of the M*-Z-SFR relation 
to be studied dir e ctly, u sing the non-parametric method 
of |Salim et al.| ( |2014| l or the related “local analog” 
method presented in this paper. Local analog method 
consists of identifying a local galaxy in SDSS (or an¬ 
other large-volume spectroscopic survey) whose stellar 
mass and SSFR match closely that of a high-redshift 
galaxy, and evaluating systematic differences between 
the metallicities of local analogs and high-redshift sam¬ 
ple (Figures and 1^. 

2. The M»-Z relation at z 2.3 shows a statistically sig¬ 
nificant dependence on SSFR at intermediate masses 
(9.7 < logM* < 10.7), the same mass range where 
such dependence (i.e., Z-SSFR anti-correlation) is seen 
in local (SDSS) galaxies with similarly high SSFRs. 
Above logM* = 10.7 no conclusions can be drawn 
because of the very small number of such galaxies 
in the sample. Anti-correlation is not present for 
lower masses (logM* ^ 9.5) KBSS galaxies, a behav¬ 
ior which appears to also hold for SDSS galaxies with 
such mass and KBSS-like SSFRs. (Figure]^. 

3. M*-Z-SFR relation of z 2.3 KBSS sample shows no 
offset with respect to the local M*-Z-SFR relation at 
lower masses (logM* < 10), regardless of whether N2 
or 03N2 lines are used to derive the metallicities (Fig¬ 
ure]^. 

4. An offset between high-redshift and local M*-Z-SFR 
relations does emerge at higher masses, reaching, 
around logM* = 10.6, ^ -0.2 dex for N2 and ~ -0.3 
dex for 03N2 metallicity. The sense of the offset is 
that z ~ 2 metallicities are lower than what is expected 
from the local M*-Z-SFR relation. M*-Z-SFR rela¬ 
tion is therefore altogether not redshift invariant (Fig¬ 
ures]^ an d|5). 


5. This high-mass offset becomes consistent for N2 and 
03N2 (-0.26 dex) if we correct some high-mass, high- 
redshift galaxies for the effects of unrecognized AGN 
contribution. AGN contamination is implicated be¬ 
cause it boosts both N2 and 03, which, when inter¬ 
preted as metallicity, leads to overestimates based on 
N2 and underestimates based on 03N2. AGN indica¬ 
tors other than the line ratios will have difficulty recog¬ 
nizing these galaxies as AGN even locally (Figures 

and|^. 

6. Current high-redshift spectroscopic surveys are biased 
against high-metallicity galaxies because they would 
have [OIII]5007 lines that fall below the [OIII] sensitiv¬ 
ity thresholds. This selection effect could in principle 
explain the gap in M*-Z-SFR relations at high mass 
and open up the possibility for a redshift-invariant rela¬ 
tion. However, the majority of current non-detections 
have weak [Nil], rather than [OIII], and few non¬ 
detections are at high mass, making this scenario for 
removing the M*-Z-SFR offset unlikely (Figures IT] |8| 
and§. 


7. Local “analogs” of high-redshift samples (selected only 
to have the same stellar mass and SSFR) are simi¬ 
larly displaced in the upper part of the BPT diagram 
with respect to the bulk of low-redshift galaxies as 
the high-redshift samples, suggesting that the lower- 
metallicity local analogs may have ISM conditions in 
common with high-redshift populations. The consis¬ 
tency between the line ratios of the high-redshift and lo¬ 
cal analogs supports our implicit assumption that equal 
line ratios at a given SSFR (and M*) indicate the same 
metallicities, even if the absolute value of this metal¬ 
licity may not be accurately given by widely-used local 
calibrations. In other words, the inapplicability of local 
calibrations to high-redshift samples is not due to the 
fact that they are local per se, but rather that they are 
largely based on typical local galaxies, and therefore 
do not account for the behavior of outlier populations 
such as the galaxies with very high SSFR for their mass 
(Figures [^and[^. 


8. The discrepancy between N2 and 03N2 metallicities 
reported at z ^ 2 (jZahid et al.|2014bl|Steidel et al.|2014 


Sanders et al.||2015| ) has three independent causes. A 

smaller part (~ 1/4) stems fr om a local mismatch in 
linear calibrations as given by |Pettini & Pagei] ( |2004| ). 
We remove th is offset by reca l ibratin g the N2 conver¬ 
sion to match IPettini & Pagel] (|2004|) 03N2 metallici¬ 
ties of local SDSS galaxies. Larger part is because of 
the changed ISM conditions at z ~ 2, as suggested in 
previous studies. However, we find that this offset is 
present to same degree in local “analogs” (Item 7). Fi¬ 
nally, the largest outliers are consistent with being ad¬ 
ditionally offset due to an AGN contribution (Item 5) 
(Figures [^and[^ . 


Additional conclusions include: 


9. KBSS 03 line ratios are on average higher than those 
of MOSDEF galaxies above logM* = 10.2, for reasons 
that may not be fully accounted for by the differences 
in the sample selection. We cannot tell which values 
should be more typical at z ~ 2 (Figure^. 
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10. Intrinsically, the AGN demarcation line at z ~ 2 proba¬ 
bly lies no more than 0.0 5 dex higher than the loc al em¬ 
pirical demarcation line ( Kauffmann et al.|2003 i. Ap¬ 
pearance of z ^ 2 BPT diagrams suggests otherwise be¬ 
cause of the scatter from larger measurement errors and 
because some of the galaxies that are not considered to 
be AGN probably are (Figure^. 


11. KBSS, and to some extent MOSDEF, have log SSFR 
versus logM* distributions (“main sequences”) that are 
steeper (more mass-dependent ) than that of z 2.3 
galaxies from 3D-HST survey ( [Whitaker et al.|[2014| l. 
In particular, at logM* >10, the 3D-HST “main” se¬ 
quence is on average ^ 0.4 dex higher than either KBSS 
or MOSDEF. While this potential bias in spectroscopic 
samples is unlikely to affect our M»-Z-SFR relation 
analysis, its sources need to be investigated (Figure [7]). 


In addition to presenting many new results, we have high¬ 
lighted observational issues that need to be fully understood in 
future work before more definitive conclusions on the subject 
of chemical evolution can be drawn. 


We thank Ryan Sanders and Alice Shapley for clarifications 
regarding their work and useful feedback on the manuscript. 
We also thank the referee Simon Lilly for constructive re¬ 
marks. 
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